Moving breathers in bent DNA with reahstic parameters 
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Recent papers have considered moving breathers (MBs) in DNA models including long range 
interaction due to the dipole moments of the hydrogen bonds. We have recalculated the value of 
the charge transfer when hydrogen bonds stretch using quantum chemical methods which takes into 
account the whole nucleoside pairs. We explore the consequences of this value on the properties 
of MBs, including the range of frequencies for which they exist and their effective masses. They 
are able to travel through bending points with fairly large curvatures provided that their kinetic 
energy is larger than a minimum energy which depends on the curvature. These energies and the 
corresponding velocities are also calculated in function of the curvature. 
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I. INTRODUCTION 

Breathers are localized oscillations in coupled networks 
of nonlinear oscillators They have been extensively 
studied in the last years 0, 0, 0, IE 0- Under some 
conditions they can move easily along the network while 
maintaining their localization Jj . A physical system 
where they may play a significant role is DNA, because 
the hydrogen bonds between nucleotides are highly non- 
linear and their openings can be related to biological pro- 
cesses such as transcription, replication and, in the case 
of MBs, with the transport of information, energy and 
charge 

DNA is also a flexible chain, and the question arises 
of to what extent MBs are hindered in their movement 
by the points of bending, or if, on the other hand, these 
points may have a biological role by trapping breathers 
and accumulating their energy. 

The shape of the DNA molecule is felt because the hy- 
drogen bond is a polar one, and the dipolar interaction 
between the dipoles is a long-range one which depends on 
the distances and orientations of the dipoles, and, there- 
fore, on the shape. 

Soliton properties in the Discrete Nonlinear 
Schrodinger Equation framework with dipole-dipole 
interaction have been considered in Rcfs. 0, IT^ . 
among others, whereas MBs in DNA Klein-Gordon mod- 
els have been studied in Ref. llSl and their interaction 
with bending points in Ref. A problem for these 



studies is the fairly large number of physical parameters 
that are not well known, among them, crucially, the 
coefficient of dipole interaction. In this paper we present 
new quantum chemical calculation that lead to the 
value of this parameter and investigate its consequences. 
Moreover, the reduction of the number of parameters 
allow us to explore the variation of another one, the 
breather frequency, which previously was considered as 
fixed. 



II. QUANTUM CHEMICAL CALCULATION OF 
THE CHARGE TRANSFER 

The only quantity that it is needed to implement the 
dipole-dipole interaction is the charge transfer q. It is de- 
fined in the following way: if po is the dipole moment of a 
Watson-Crick hydrogen bond at the equilibrium distance, 
and it is stretched by a small amount u, the new dipole 
moment is given to the first order in u hy p = pq + qu . 
The relevant estimations have been obtained in Ref. |l3j | , 
yielding q values between — 0.0014 e and — 0.0183 e for 
an A-T base pair and between — 0.025 e and — 0.055 efor 
a G-C base pair, depending on the quantum-chemical 
method used. However, these values were obtained only 
for Watson-Crick base pairs in vacuo. In the present 
work we significantly extend our model, in that we take 
into account 

a) the variation of the dipole moment of the whole 
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nucleoside (base + deoxyribose) pair, 

b) the influence of the DNA duplex environment on 
dipole moments of the nucleoside pairs in question. 

We also consider the helical structure of the DNA 
double strand. The procedure now is as follows: we 
consider regular homogeneous stacked trimers of nucle- 
oside pairs, namely, adeninosine-thymidine (AT) and 
guanosine-cytidine (GC) pairs in the following arrange- 
ment: AT/AT/AT and GC/GC/GC. We take the 'effec- 
tive' dipole moment of a nucleoside base pair to be one 
third of the trimer dipole moment. Then we stretch and 
squeeze H-bonds by gradually adding or subtracting of 
up to 0.1 Angstrom to the equilibrium spacing between 
A and T (or G and C) in the central pair of the trimer, 
with the two flanking pairs remaining in the standard 
B-DNA-conformation. At each of the H-bond stretch- 
ing/squeezing states, the 'effective' dipole moment of a 
nucleoside pair was estimated using semiempirical quan- 
tum chemistry. The relevant computational details are 
published elsewhere 0|. 

We evaluate a linear regression of the 'effective' dipole 
moment (p) onto the nucleotide spacing change (u) in the 
Watson-Crick pairs {p — po + qu), where the coefficient 
q, delivers the desired estimate. As a result, we have 
obtained po = 1.84 D, q = — 0.09e, for the AT pair and 
Po = 3.15 D, q = -0.10 e, for the GC pair. The Watson- 
Crick H-bond stretching diminishes the effective dipole 
moment of nucleoside pairs, which ought to be connected 
with the proper changes in the purine-pyrimidine molec- 
ular orbital overlap included into quantum-chemical eval- 
uations of dipole moments. If we take into account only 
the conventional electrostatic charge distributions, in- 
creasing the spacing between the partners in the Watson- 
Crick base pair should lead to the increase in the dipole 
moment. Note that the values of q obtained in this 
way are substantially larger than the ones obtained in 
Ref. 13] and very similar for both the A-T and G-C 
pairs, in spite of the different number of hydrogen bonds. 
This could be explained by the appreciable influence of 
the DNA duplex surrounding on the Watson-Crick A- 
T and G-C base pair dipole moments. We agree that 
taking an "effective" dipole moment of the stacked base 
pair trimer is only an approximation. We intended to 
include the coupling between the nearest neighbours in 
the DNA stack. That the "effective" dipole moment val- 
ues differ from those of the isolated base pairs is actu- 
ally a signiflcant result showing that we were right to 
include the DNA stacking interactions at least approx- 
imately into our model. As to the Coulomb electron- 
electron interaction, it was considered at our Hartree- 
Fock-level quantum-chemical calculations in the form of 
Pauli correlations, where electrons of the same spin are 
repellent. Note also that being the values of q for both 
homopolynucleotide duplexes very similar we should not 
expect q to change significatively for heterogeneous DNA. 



III. THE MODEL 

The Hamiltonian system also used in Ref. 0, Q| is 
a Peyrard-Bishop model 16] augmented with long-range 
interaction. It can be written as: 



1 

n=l 

(w„+i - -u„)2 + ^ X! ■ 

i 

The explanation of the terms and variables is as fol- 
lows: Un is the stretching of the n-th hydrogen bond; 
m is the reduced mass of a nucleotide pair; the Morse 
potential D (e~^"" — 1)^, represents the energy of a hy- 
drogen bond between a pair of nucleotides; (m„+i — 
is the stacking energy between neighboring bonds; 
i X^i JinUi Un is the energy of the dipole-dipole interac- 
tion; if r„ is a vector that denotes the position of each 
particle in the 2-d space, Jin = J/\rn — ri\'^, for i ^ n 
and otherwise, where J, hereafter referred to as the 
dipole parameter, is given by J = g^/(4 Treo c?'^), q and 
d — 3.4 A, being the charge transfer and the distance be- 
tween neighbouring base pairs, respectively, as deduced 
in Ref. jlj]- We approximate the shape of the molecule 
in the vicinity of a bending point by a parabola with cur- 
vature K embedded in a plane, and consider the dipoles 
with the same orientation and orthogonal to it. We con- 
sider fixed geometry as the H-bond vibrations are much 
faster than the bending movement. Thus, r„ = {xmVn) 
with Un ~ K.x\l2. We also neglect the heterogeneity of 
DNA. 

The values of the parameters are taken from Ref. ^3 ■, 
i.e., D = 0.04 eV, h = 4.45 A'^ and to = 300 amu. 
As discussed in the same reference, fc can take values be- 
tween 0.01 eV/A2 and 10 eV/A^ and, m consequence, we 
consider it as an adjustable parameter in that range. The 
value of J for a charge transfer q intermediate between 
the A-T and G-C base pairs is J = 0.0031 eV/A^. 



IV. MOVING BREATHERS IN A STRAIGHT 
CHAIN 

Breathers, exact to machine precision, are calcu- 
lated numerically using techniques based in the anti- 
continuous limit p^ . MBs with good movability prop- 
erties, i.e., low dispersion, are obtained for values of the 
parameters in the vicinity of a bifurcation point, called 
inversion of stability, where a single breather becomes 
unstable and a double one stable, or vice versa, and per- 
turbing them with a vector coUinear to the eigenvector 
that becomes unstable 0, H] ■ 

The particles of the breather oscillate coherently with 
a frequency v^. For the sake of simplicity, we define the 
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FIG. 1: Range of existence of moving breathers. The upper 
limit corresponds to the bifurcation of stability inversion. The 
lower limit is the top frequency of the phonon band. The 
dispersion is smaller if closer to the upper limit. 
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FIG. 3: Dependence of the critical translational energy Kc 
with respect to k for k = ks- The frequencies fib of the 
breathers are (from up to down) 0.95, 0.90, 0.85, and 0.80. 
Smaller values of ft are not shown as Kc is very small and no 
accurate results can be obtained. 




FIG. 2: Dependence of the effective mass m* with respect to 

the breather frequency for three different values of the cou- FIG. 4: Same as Fig.Olbut for the critical initial translational 
pling constant k: 0.2376, 0.2693 and 0.3010 eV/A^ from bot- velocity (vc). 
torn to top. 



non-dimensional frequency fib = Vh/vo, where i/q, de- 
fined through 



'^o^T^V = T13THz (2) 

ZTT V m 

is the linear frequency of the isolated oscillators in Eq. ^ 
As the on-site potential is soft, fib < 1- We have lim- 
ited its values to the interval (0.65, 0.95) because lower 
frequencies brings about breathers difficult to move, due 



to their small width, and MB with larger ones develop 
frequencies which interact with the phonon band. The 
latter is due to the fact that MBs are a wave packet 
with frequencies around the one of the perturbed static 
breather. 

For a given value of the breather frequency Sib, MBs 
exist for values of k above a critical one kg = fcs(^'b), for 
which the inversion of stability takes place. This value 
has to be calculated numerically. 

The phonon band is composed of the frequencies of the 
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linear modes, given by 0, 0| : 

where q is the wave number, lUq — 2 7: and Fs{z) — 
"Y^j^yi jk^ is the Polylogarithmic or Jonquiere function. 
They exist for increasing values of fc > fcg until an upper 
limit km = km{vh), for which the top of the phonon band, 
at g = TT, reaches the second harmonic of the breather. 
km is given by 

km^(^l--^mLul^^jai): (4) 

where Q{s) is the Riemann's zeta-function. This is also 
the upper limit for MBs existence. Fig.nshows the range 
of existence of MBs. It can be seen that they exist for 
a wide range of frequencies, but they have much smaller 
dispersion at the vicinity of the inversion stability curve 
and we will restrict often our study to them. 

Static breathers in the above mentioned region have 
an energy between KbT and 3.7 KbT, with T = 310 K. 
The dipolar energy of these breathers oscillates between 
0.006 KbT and 0.023 KbT. The energy depends mono- 
tonically with the frequency, corresponding the highest 
energy to the smallest frequency. 

Moving breathers behave as quasiparticles, having an 
effective mass m* related with the kinetic energy of the 
breather by if = \Tn* v"^ , with v being the translational 
velocity of the breather 0,13. K e (1.5, 50) x 10-* KbT is 
also the kinetic energy added to the static breather, m* is 
approximately constant as long as v is small enough and, 
therefore, emits small phonon radiation. Fig.Elshows the 
dependence of the mass with respect to the frequency for 
three different values of k. 

Thus, MBs have an effective mass m* between 30 and 
46 amu, depending on the value of frecuency and the 
stacking parameter. 

V. MOVING BREATHERS IN BENT CHAINS 

In order to study the effect of the curvature on breather 
mobility, we launch a breather through the ben ding 
point. The bending acts as a potential barrier [lj| 
in accordance to the trapping hypothesis formulated in 



Ref. ^3 (fo'" ^ discussion of this last point see Ref. 01 )■ 
It implies that breathers are reflected when reaching the 
bending point as long as their translational energy K is 
below a critical point Kc- Fig.Oshows the dependence of 
this value with respect to k for several frequencies. The 
dependence of the critical initial velocity is also shown 
in Fig. 21 Kc, which is also the amount of energy given 
to a static breather to start its movement, is fairly small 
compared to the thermal energy at physiological tem- 
perature of 0.027 eV. Therefore, it seems likely that it 
can be given easily by the environment, and most MBs 
are able to pass through bending points without being 
reflected. 

VI. CONCLUSION 

In this paper we consider a model for DNA with long- 
range interaction due to the dipole moments of the hy- 
drogen bonds among nucleotides. The dipole parameter 
J has been obtained through quantum chemical calcula- 
tions including the whole nucleoside. The consequences 
for moving breathers in straight and bent DNA chains 
have been analyzed. These can be summarize as: 1) For 
each different value of the stacking parameter there ex- 
ist MBs with different frequencies, including an upper 
limit, determined numerically for the inversion of stabil- 
ity curve, and a lower limit determined both numerically 
and analytically by the top of the phonon band. 2) MBs 
breather have a fairly large effective mass of 30-50 amu. 
2) The values of J, though small, are large enough for 
the MBs to feel the the bending of DNA. 3) We have 
calculated the minimum kinetic energy and the velocity 
of a breather to be able to pass over a bending point in 
function of the curvature. As the minimum kinetic en- 
ergy is small, it seems plausible that most MBs will have 
enough energy to pass through bending points, although 
the ones with low energy will be reflected. 
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